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Abstract 

In this paper we study the one zero textures of low energy neutrino mass matrices in presence of 
a sterile neutrino. We find that the mass matrix elements m a p (a, (3 = e, fi, r) involving only the 
active states can assume vanishing values in the allowed parameter space for all the mass spectrum. 
Among the mass matrix elements connecting the active and sterile states, m es and m^ s can become 
small only for the quasi-degenerate neutrinos. The element m rs on the other hand can be very 
small even for lower values of masses since the 3-4 mixing angle only has an upper bound from 
current data. The mass matrix element (m ss ) involving only the sterile state stays ~ 0(1) in the 
whole parameter region. We study the possible correlations between the sterile mixing angles and 
the Majorana phases to give a zero element in the mass matrix. 



" Email Address: monojit@prl.res.in 



'fEmail Address: sruba@prl.res.in 

■^Email Address: shivani@cski mTyonsei.ac.kr| 



3 Email Address: cskim@yonsci.ac.kr 



1 



I. INTRODUCTION 



Light sterile neutrinos were invoked to explain the results of the LSND experiment which 
reported escalation events in the », - „ e mode corresponding to a mass squared difference 
~ eV 2 [l|. Adding one sterile neutrino to the standard 3 generation framework gives rise to 
two possible mass spectra - the 2+2 in which two pairs of mass states are separated by a 
difference ~ eV 2 and 3+1 in which a single predominantly sterile state differs by ~ eV 2 from 
the three active states ^|. Subsequently the 2+2 schemes were found to be incompatible 



with the solar and atmospheric neutrino data |3J. The MiniBoone experiment was designed 

to test this and its antineutrino data confirmed the LSND anomaly [4J. Both 3+1 and 3+2 

sterile neutrino schemes have been considered to explain these results js — 7 1 - Such global 

fits aim to explain the non-observance of eV 2 oscillations in the disappearance channel in 

other short baseline experiments as well as the reported evidence in LSND/ MiniBoone 

experiments. The relevant probabilities for 3+1 case is governed by a single mass squared 

difference and hence is independent of the CP phase. In the 3+2 scheme dependence on CP 

phase comes into play and one gets a slightly better fit. 

Other evidences in support of sterile neutrinos include - the reactor and Ga anomaly. The 

first one refers to the deficit in the measured electron antineutrino flux in several experiments 

when the theoretical predictions of reactor neutrino fluxes were reevaluated [8|. The second 

one implies shortfall of electron neutrinos observed in the solar neutrino detectors GALLEX 

f] 

and SAGE using radioactive sources [9[. Both these can be explained by adding light sub-eV 
sterile neutrinos in the three generation framework. 

There has also been some hint in favour of sterile neutrinos from cosmological observations 
of a " dark radiation" which is weakly interacting and relativistic. Attributing this to sterile 
neutrinos one gets the bound on the number of neutrinos as N e ff = 4.34 ± 0.87 at 68% C.L 



10| . The Plank satellite experiment which has very recently declared its first results |ll| . 



on the other hand, gives N e ff = 3.30 ± 0.27 at 68% C.L. which allows for an extra sterile 
neutrino at 95% C.L., although its mixing with active species can be very tightly constrained 



121 ] within the framework of standard cosmology. Thus the sterile neutrinos continue be 



intriguing and many new experiments are planned proposed to test this [13]. 

Theoretically, sterile neutrinos are naturally included in Type-I seesaw model fl^ . But 
their mass scale is usually very high to account for the small mass of the neutrinos. Light 
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sub-eV sterile neutrinos as suggested by the data can arise in many models |l3j ]. 

Irrespective of the mechanism for generation of neutrino masses the low energy Majorana 
mass matrix in presence of an extra sterile neutrino will be of dimension 4x4 with ten 
independent entries and is given as, 

M v = V*M^ a9 V ] (1) 

where, M^ mfl = Diag(mi,rri2, 7713,7714) and V denotes the leptonic mixing matrix in a basis 
where the charged lepton mass matrix is diagonal. One of the important aspects in the 
study of neutrino physics is to explore the structure of the neutrino mass matrices. At 
the fundamental level these are governed by Yukawa couplings which are essentially free 
parameters in most models. These motivated the study of texture zeros which means one 
or more elements are relatively small compared to the others. Texture zeros in the low 
energy mass matrices in the context of three generations have been extensively explored 



both in the quark and lepton sector [15J], [16|. Such studies help in understanding the 



underlying parameter space and the nature of the mass spectrum involved and often predict 
correlations between various parameters which can be experimentally tested. For three 
generation scenario it is well known that the number of maximum texture zeros in low 
energy mass matrix is two [15 1. In the context of the 4-neutrino case however more than two 

n 

zeros can be allowed [171 ] . Two zero textures of sterile neutrinos have been studied recently 
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171 ] and three zero cases have been considered in 181 ] . In this paper we concentrate on 
the textures where one of the mass matrix elements is vanishing. For the 4x4 symmetric 
mass matrix it gives total 10 different cases which needs to be investigated. We study the 
implications of one zero textures and the possible correlations between the parameters. We 
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20]. 



also compare our results with the 1 zero textures for three active neutrinos 

The plan of the paper goes as follows. In section II we discuss the possible mass spectra 
and the mixing matrix in the 3+1 scenario. In the next section we present our study 
regarding the implications of one vanishing entry in the low energy neutrino mass matrix. 
We conclude in section IV. 
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Figure 1: The allowed 3+1 mass ordering 



II. MASSES AND MIXING IN THE 3+1 SCHEME 

There are two ways in which one can add a predominantly sterile state separated by ~ 
eV 2 from the standard 3 neutrino mass states. In the first case the additional neutrino can 
be of higher mass than the other three while in the second case the the fourth neutrino is the 
lightest state. The later turns out to be incompatible with cosmology since in this case three 
active neutrinos, each with mass ~ eV results in an enhanced cosmological energy density. 
Thus it suffices to consider only the first case which admits two possibilities displayed in 
Fig. 1. 

(i) SNH: in this rrii rs m 2 < m 3 < m 4 corresponding to a normal hierarchy (NH) among 
the active neutrinos which implies, 

m 2 = \f m\ + Amf 2 , = a/ mf + Amf 2 + Am^ , = \/m\ + Amf 4 . 

(ii) SIH : this corresponds to m-^ < ni2 ~ vn\ < implying an inverted ordering among 
the active neutrinos with masses expressed as, 

mi = a/ m| + Am 2 3 , m 2 = \/ m| + Am^ + Am 2 2 , 7714 = ^/m^ + Am| 4 . 



4 



Here, Am| = m| — mf. In the extreme cases these masses can be written as 



SNH : |m 4 | 



Am^ > |m 3 | 



Am^ 



3> m 2 



Am| 2 ^> | 



(2) 



S7iJ : |m 4 | 



Ami > |m 2 | 



mi 



Am^ > |m 3 | 



(3) 



SQD : |m 4 | ^> |mi| 



m 2 



m . 



(4) 



The first two cases correspond to complete hierarchy among the active neutrinos while 
the last one is the quasi-degenerate (QD) regime where the three active neutrinos have 
approximately equal masses. 

In the 3+1 scenario, the neutrino mixing matrix, V in the flavor basis can be parametrized 
in terms of sixteen parameters. In addition to the three mixing angles between the active 
flavors, (#i3, #i2, #23) we now have three more mixing angles from sterile and active mixing, 
(#14, #24, #34)- There are six CP violating phases, three Dirac (<5i 3 , 5 U , 5 U ) and three 
additional Majorana phases as (a, (3, 7) as neutrinos here are considered to be Majorana 
particles. Then, there are four masses of neutrino mi, m 2 , m 3 corresponding to three active 
states and which is predominantly the mass of heavy sterile neutrino. 

The mixing matrix V can be expressed as V = U.P 2l| where 



U — R34R24R14R23R13R12 



(5) 



where Rij denotes rotation matrices in the ij generation space and is expressed as, 
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Here we use the abbreviations Sjj = sin #jj and Cjj = cos#jj. The phase matrix is diagonal 
and is expressed as, 



P = Diag(l, e ia , e i(/3+5l3) , e i(7+<5l4) ). 
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Parameter 
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3cr range 


Amf 2 /10" 5 eV 2 (NH or 
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7.54 


6.99 - 8.18 


sin 2 (912/KT 1 (NH or IH 


) 


3.07 


2.59 - 3.59 


Am|o/10- 3 eV 2 (NH) 




2.43 


2.19 - 2.62 


Amf^/lO" 3 eV 2 (IH) 




2.42 


2.17 - 2.61 


sin 2 e n /m- 2 (NH) 




2.41 


1.69 - 3.13 


sin 2 tfis/lO" 2 (IH) 




2.44 


1.71 - 3.15 


sin 2 ^/lO^ 1 (NH) 




3.86 


3.31 - 6.37 


sin^a/KT 1 (IH) 




3.92 


3.35 - 6.63 


Am? SJV n ( Am? 4 or Ami 


) eV 2 


1.62 


0.7 - 2.5 


sin 2 6*i4 




0.03 


0.01 - 0.06 


sin 2 024 




0.01 


0.002 - 0.04 


sin 2 i?34 






< 0.18 


C/10" 2 (NH) 






2.7 - 3.7 


C/1(T 2 (IH) 






2.7 - 3.8 


£/10 3 (NH) 






0.27-1.14 


£/10 3 (IH) 






0.27- 1.15 



TABLE I: 3<r ranges of neutrino oscillation parameters 
neutrino parameters are from 
and £ 



23] 



221 ] . The current constraints on sterile 



241 ] . Also given are the 3a ranges of the mass ratios £ 



Amf 2 
|Am| 3 | 



Am lsND 



The best-fit values and the 3a ranges of the oscillation parameters in the 3+1 scenario are 
given in Table I where in addition to the masses and mixing angles we also present the mass 

AtTI^ /\m ^ 

ratios £ = t-t — ti and t = . LS i'P which would be useful in our analysis. 

Z\m 23 | |Am 23 | 

In Fig. 2 we have plotted the sum of neutrino masses against the lowest neutrino mass 
for both NH and IH. The band corresponds to variation of the mass squared differences in 
their current 3a range. We also show the cosmological upper bound on neutrino masses in 
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Figure 2: Plots of sum of light neutrino masses (Em„) vs the lowest mass in the 3+1 scenario. 
Also shown are the the cosmological upper bound on neutrino mass from the analysis of CMB 
data plus matter power spectrum information (SDSS) and a prior on //o(HST) and from the 
analysis of CMB data only from reference [23j] for 3+1 scheme. 



3+1 scenario from 23j. The combined analysis of CMB + SDSS + HST seems to rule out 
the mass spectrum of 3+1 scenario whereas if only CMB data is taken then region for the 



lowest mass < 0.4 eV gets allowed for both the hierarchies. Note that the analysis in [23 1 
does not incorporate the Planck results hJ which can constrain the sum of masses further, 
In our analysis we have varied the lowest mass up to 0.5 eV. 



III. NEUTRINO MASS MATRIX ELEMENTS 



In this section we study the implication of the condition of vanishing m a p for the 3+1 
scenario, where a, (3 = e, /x, r, s. Since m a/ j is complex the above condition implies both real 
and imaginary parts are zero. Therefore to study the 1-zero textures we consider \m a p\ = 0. 

A. The Mass Matrix element m ee 

The matrix element m ee in the 3+1 scenario is given as, 

2 2 2 2 2 2 2ia , 2 2 2iB , 2 2ri fr\ 

mee — m\c 14 c 13 c 12 + m 2 s 12 c 14 c 13 e + m 3 s 13 c 14 e p + m 4 s 14 e y . (6) 
This is of the form 

m ee = c 2 u (m ee ) 3l/ + e 2i7 Si 4 m 4 , (7) 
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Figure 3: Plot of \m ee \ versus the lowest mass. The panel (a) corresponds to the three generation 
case while panel (b) is for 3+1 case. In panel (b) all the mixing angles are varied in their 3<t 
range and the CP violating phases are varied in their full range. The panel (c) and (d) are for 
specific values of #14 with all other parameters covering their full range. 

where (m ee )s u corresponds to the matrix element in the 3 active neutrino case. The contribu- 
tion of the sterile neutrino to the element m ee depends on the mass 7714 and the active-sterile 
mixing angle #14. Of all the mass matrix element m ee has the simplest form because of the 
chosen parametrization and can be understood quite well. Using approximation in Eq. ([2]) 
for the case of extreme hierarchy one can write this for NH as, 



m ee » c 2 1A (m ee ) 3u + e 2ll s 2 u J Am?. 



where (m, 



ee)3u 



Am 



23 V 



3 2ict„2 „2 
L 13 6 12 



2 2ifi 
*13 e 



'14i 



Am? 



(8) 



and, C = iX" i 2 \ is t ne ratio of the solar 

|Am 23 | 

and the atmospheric mass squared differences. The modulus of m ee is the effective mass 
that can be extracted from half life measurements in neutrinoless double beta decay. In 
Fig. 3 we plot the effective mass as a function of the smallest mass by varying #14 in its 
complete 3<r range from Table I as well as for specific values of the mixing angle #14. The 
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Majorana phases are varied randomly in the range to tt in all the plots. The first panel is 
for #i4 = i.e the three generation case. It is seen that for present values of the oscillation 
parameters the cancellation condition is not satisfied for mj — > for NH. However, as one 
increases m 1; complete cancellation can be achieved. For IH the complete cancellation is 
never possible. These results change when we include the sterile contribution as is evident 
from the panel (b) in Fig. 3 which shows the effective mass for NH and IH by varying all 
the parameters in their full 3a allowed range. The behaviour can be understood from the 
expressions of \m ee \ in various limiting cases. For NH, in the hierarchical limit of mi — > 
the major contributor will be the additional term due to the sterile neutrinos because of 
higher value of m 4 . Complete cancellation is only possible for smaller values of #14 so that 
this contribution is suppressed. The typical value of #14 required for cancellation can be 
obtained by putting a = (3 = (which would maximize the three neutrino contribution) 
and 7 = 7r/2, as 

tan 2 6u w (vfo; 3 s|+sj 3 ) w 10 _ 3) (9) 



which lies outside the allowed range of 6*14 given in Table I. Here f = tt — rr. As we increase 
777,1, (m ee ) 3l , increases and can be of the same order of magnitude of the sterile term. Hence 
one can get cancellation regions. The cancellation is mainly controlled by the value of 8 U . 
For higher values of s 2 4 one needs a higher value of m x for cancellation to occur. The 
correlation between ni\ and #14 is brought out by the panels (c) and (d) in Fig. 3. 

For IH case, in the limit of vanishing m 3 using approximation in Eq. ([3]), m ee in a 3+1 
scenario can be written as 



(:'(:' I 



144^4(4 + s 2 12 e 2 n + ^4,sle^\. (10) 



The maximum value of this is achieved for a = 7 = which is slightly lower than that of 
NH in this limit. The element vanishes in the limit m 3 w eV when a = and 7 = 7r/2 
provided 



c 2 

tan 2 6 U w w 0.05 (11) 
This is well within the allowed range. This behaviour is in stark contrast to that in the 3 



neutrino case 25[ . There is no significant change in this behaviour as the smallest mass 
m 3 is increased since this contribution is suppressed by the s 2 3 term and the dominant 
contribution to (m ee ) 3u comes from the first two terms in Eq. ([6]). Therefore in this case we 



do not observe any correlation between m 3 and s 2 



ii- 



While moving towards the quasi-degenerate regime of mi m 777 2 ~ WI3 we find that 
effective mass can still be zero. However, when the lightest mass approaches a larger value 

0.3 eV we need very large values of active sterile mixing angle #14, outside the allowed 
range, for cancellation. Hence the effective mass cannot vanish for such values of masses. 

Also shown is the current limit on effective mass from combined KamLAND-Zen and EXO 
200 results on the half-life of 0^/3/3 in 136 Xe Q, 0. When translated in terms of effective 
mass this corresponds to the bound \m ee \ < 0.11 — 0.24 eV including nuclear matrix element 
uncertainties. For the three generation case, the hierarchical neutrinos cannot saturate this 
bound. But in the 3+1 scenario this bound can be reached even for very small values of 777.3 
for IH and for some parameter values it can even exceed the current limit. Thus from the 
present limits on neutrinoless double beta decay searches a part of the parameter space for 
smaller values of 7773 can be disfavoured for IH. For NH, the KamLAND-Zen + EXO 200 
combined bound is reached for mi = 0.02 eV and again some part of the parameter space 
can be disfavoured by this bound. 

B. The Mass Matrix element m efl 

The mass matrix element m efl in the presence of extra sterile neutrino is given as 

m eil = c^e^-^+^Wu^ + e i(<5l3+2/3 W 13 (c 13 c 2 4S 2 3 ~ e^ 14 "**-*^ (12) 

s 13 s uS2a) + c 12Cl3777 1 (— C23C 24 Si2 + C 12 (— e %&1A C24S13S23 — e^ 14 ^ C13S14S24)) 

+ e 2ia 777 2 ci3Si2(ci2C 2 3C24 + s 12 (~e iSvi C24S13S23 - e i(5l4 ~ 524) ci 3 si 4 S24) ) ) • 

Unlike 777 ee here the expression is complicated and an analytic understanding is difficult from 
the full expression. The expression for m eil in the limit of vanishing active sterile mixing 
angle # 2 4 becomes 

«7 eA1 = ci 4 (m eAt ) 3l ,. 

Since the active sterile mixing is small, in order to simplify these expressions we introduce 
a quantity A =0.2 and define these small angles to be of the form aX. Thus a systematic 
expansion in terms of A can be done. For sterile mixing angle 

sin0i4 « 6> 14 = xuA, 
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sin6>24 ~ 6*24 = X24A, (13) 

and the reactor mixing angle as 

sin 6»i3 « 6»i3 = X13A. (14) 

Here Xij are parameters of 0(1) and their 3a range from the current constraint on the mixing 
angles is given by 

Xis = 0.65-0.9, (15) 
Xu = 0.5-1.2, 
X24 = 0.25 - 1. 

Using the approximation in Eqs. ([2]), (1131) and (1T41) we get the expression for |m eM | for 
normal hierarchy as 



'Am 2 3 {y/(s 12 c 12 c 23 e 2ia + e^(e 2i ? - e 2ia y/Cs 2 2 )s 23 X X i3 (16) 
+ A 2 e i(^-^) (e 2n^_ e ^^4) XuX24 }|. 

To see the order of magnitude of the different terms we choose vanishing Majorana phases 
while Dirac CP phases are taken as it. The mass matrix element m eiI vanishes when 

^■512012023 - (1 - V / C S 12) S 23AXl3 + A 2 ( ~ V^OXhXM = 0. (17) 

The three generation limit is recovered for s\ 4 = and in panel (a) of Fig. 4 we show \m e/Ji \ 
as a function of m\ of this case, for NH. Panel (b) (red/light region) of Fig. 4 shows |m eM | 
for the 3+1 case, with all parameters varied randomly within their 3a range. The figures 
show that \m efl \ = can be achieved over the whole range of the smallest mass for both 3 
and 3+1 cases. However, we find that in the hierarchical limit cancellation is not achieved 
for large values of #24, since in that case the third term of Eq. (HTJ) will be of the O (10" 1 ) 
compared to the leading order term which is of the O (10~ 2 ) and hence there will be no 
cancellation of these terms. This can be seen from panel (b) (green/dark region) of Fig. 4 
for s 2 4 = 0.04. In the QD limit the contribution from the active terms are large enough to 
cancel the sterile contribution and thus \m efJ \ = can be achieved. 

For IH using the approximation Eq. ([3]) for the hierarchical limit we get the expression 
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Am 2 3 {c 12 s 12 c 23 (e 2ia - 1) - e^{c\ 2 + s 2 12 e 2ia )s 23X i 3 \ (18) 
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(a) (b) 

i h , — , — . — — ^ 1 h , — 




Lowest Mass(eV) Lowest Mass(eV) 

Figure 4: Plots of | function of the lowest mass mi for NH. Panel (a) correspond to the 

three generation case while (b) (red/light region) is for 3+1 case and also for s^ A = 0.04 
(green/dark region). All the parameters are varied in their full 3a allowed range and the CP 
violating phases are varied from to ir unless otherwise stated. 

- e^- 5 ^\\ lAX 2,{cl 2 - e*>y/£ + e 2ia s 2 12 )}\. 

To see the order of magnitude of the various terms we consider the case when Majorana 
phases vanish and the Dirac phases assume the value ir. Then we get for vanishing m eM , 

S23AX13 - A 2 (l - v^)Xi4X24 = 0. (19) 

In panel (a) of Fig. 5 we display the plot of \m efi \ with m 3 for the 3 generation scenario 
i.e for sin 2 6*24 = for IH. In panel (b) we consider the 3+1 case with all the parameters 
varying in their allowed range. Note that in the small limit (cf. Eq. 3.13) for a = the 
leading order term vanishes. For this case, for large active sterile mixing angle #24, the A 2 
term becomes large O (1CT 1 ) and the cancellation with A term is not be possible. When CP 
violating phase a is non zero, the leading order term can cancel the A 2 term even for large 
values of s 2 4 . These features are reflected in panel (c) where we plot |m e/ J for s| 4 = 0.04 and 
a = (blue/dark region) and by varying a in its full range (cyan/light region). As expected, 
for a = 0, cancellation is not achieved for smaller values of m^. Thus the condition \m eil \ = 
implies some correlation between m 3 and a for IH. Even if a is varied in its full range, the 
absolute value of the matrix element |m e/i | can vanish only if the product X14X24 is small, i.e. 
s 2 4 and S24 are simultaneously small . This is because if they are large the A 2 term becomes 
of the O (10" 1 ) and hence cancellation will not be possible. This is seen from panel (d) 
where for s 2 4 = 0.06 and s 2 4 = 0.04 the region where m 3 is small gets disallowed. Taking 
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Figure 5: Plots of \m efl \ vs m.3 for inverted hierarchy for (a) three generation case (b) 3+1 case 
with all parameters varied randomly in their full range. Panel (c), (d) and (e) are for specific 
values of and a. The panel (f) shows the correlation between a and when all other 

parameters are randomly varied. 

CP violating phase a = ir/2 makes the magnitude of leading order term (si2Ci2C23\/C) quite 
large and smaller values of #24 cannot give cancellation even for large values of m 3 which 
can be seen from panel (e) of Fig. 5. For the occurrence of cancellation s\ A has to be > 0.01 
for a = 7r/2 as can be seen from panel (f) where we have plotted the correlation between a 
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and sl A for \m eu \ = 0. 



The Mass Matrix element m e 



The mass matrix element m er in the presence of an extra sterile neutrino is 

m eT = c 14 c 24 e l(27+5l4 Wi4S34 + m 3 c u srse i ^ +5 ^(-c 24 s 13 s 14 s 3i e i ^ 5 ^ 

+ ^3(023034 — e i<524 S23S24-S34)) + ^n2S\ 2 Ci 3 Ci4e 2la {Ci 2 { — C 3 4S 23 — C23S24<S34e l<524 ) 

+ s 12 (-c 13 c 2 4S U s 3A e i5li - e iSl3 s 13 (c 23 c 34 - e iS24 s 23 s 24 s u ))) 

+ miCi2Ci3C 14 (-si2(-C34S23 - c 23 s 2A s 3A e %& ' 2i ) + c 12 (-ci 3 c 24 sus u e lSl4 - e lSl3 si 3 

(C23 c 34 — e* 524 S23S24S34)))- 



(20) 



This element m eT and m efl are related by \x — r permutation symmetry 

( I o\ 

10 

10 

\0 1 / 

in such a way that 



P, 



IMT 



pT • p 



For three active neutrino case the mixing angle #23 i n the partner textures linked by fi — t 
symmetry are related as /9 2 3 = (f — #23). However, in the 3+1 case the relation of #23 between 
two textures related by this symmetry is not simple. The active sterile mixing angles #24 
and 6*34 are also different in the textures connected by ji — r symmetry and are related as 
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6-10 = 8 



12, 



#13 = 9 



13, 



9-1 4 = 



14, 



sin 9 



sin #24 = sin 6*34 cos # 2 4, 
cos 6*23 cos #34 — sin 6 23 sin 6 34 sin 6 2i 



23 



cos 61a sin 6>f 4 



sin 9 



'24 1 

sin 6 2 a 



31 



1 — cos #24 sin 6>f 4 



(21) 
(22) 
(23) 

(24) 



Due to these relations the behaviour of m efl is different from that of m eT unlike in three 
active neutrino case where the plots of these two elements were same except for #23 which 
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differed in octant for both the textures. 

It is found that in the limit of small 6*24 the two active sterile mixing angles #24 ~ #34 from 
Eq (I22p . The same can be seen from Eq (1241) which gives #34 ~ #24 for smaller values of the 
mixing angle #34. Thus, for these cases the behaviour shown by #24 m i m fifi) is same as 
shown by #34 in m er (m rr ). 

In the limit of vanishing active sterile mixing angle #34 this element becomes 

m eT = c u (m eT ) 3u . 

Using the approximation in Eq (TS]) for NH the above element can be expressed as, 

\m eT \ « \yj ' Am 2 23 {-s 12 s 23 c 12 c u ^e 2ia + \(c 23 c u e t{2 ^ +Sri) X13 - c 23 c u s 2 12 (25) 
Xl3 yCe i(2Q+<5l3) + e i(2 ^ +5l4) ^34Xi4 - e^^v^Wi* ~ c 12 c 23 
e^^s 12 s uX 2,VC) ~ e^+^(e^ - e 2 ^L ^23*34X13X24^}!. 

For a case of vanishing Majorana phases and Dirac phases having the value n this element 
can vanish when 

- Ci2C 34 a/Csi2S23 ~ (1 ~ V^L^^A^IS^ + A(-C 23 C34Xl3 + (26) 
V / C C 23Sl2(c 3 4Sl2Xl3 + C12S34X24) + Sl2 S 34Xl4 V'C ~ V^^Xu) = 0. 

For a vanishing active sterile mixing angle #34 one recovers the 3 generation case. In this 
limit, from Eq. (l26p one observes that the leading order term and the term with A are of 
the same order ~ O (10~ 2 ) while the A 2 term vanishes and hence cancellation is possible 
excepting for very low values of the lightest mass. We can see this in panel (a) of Fig. 
6. In panel (b) (red/light region) all the parameters are varied randomly (3+1 case) and 
cancellation is seen to be possible over the whole range of mi. In panel (b) (green/dark 
region) we also plot the element \m eT \ for the upper limit of s 2 4 = 0.18. In this case there is 
no cancellation for very low values of the smallest mass. This is because when s\ A is large, 
the A term containing £ becomes large O (1) and there will be no cancellation. 
For inverted hierarchy the element m eT using the approximation in Eq ([3]) becomes 

\m eT \ « \^A^ 3 {c 12 c u s 12 s 23 (-e 2ta + 1) + e^ +s ^(c 2 12 + e 2 "^)^^ 2 ^^ 

- A(c23C34Xi3e^ 3 (c 2 2 + e 2ia s\ 2 ) + e iS ™ s 3 4Xi4(c 2 2 + e ia s\ 2 ) 

- e^' + ^h MXu V^ + Ci2C 23 s 12 s M X2Ae l524 (e 2 ^ - 1))}| (27) 
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Figure 6: Plots of \m er \ for normal hierarchy with lowest mass m 4 . The panel (a) corresponds to 
three generation case. In (b) (red/light region) all the parameters are varied in their full allowed 
range and the green/dark region is for s| 4 = 0.18 with all the other parameters covering their full 

range. 

In the limit of vanishing Majorana phases and Dirac CP violating phases equal to it this 
element becomes negligible when 

A(c 2 3C 34 Xl3 + <S 34X14 ~~ s 34Xl4 V£) + S23S34X13X24A 2 = 0. (28) 

In panel (a) of Fig. 7 the three generation case is reproduced by putting s| 4 = and in (b) 
all the parameters are varied in their allowed range (3+1 case). In both the figures we can 
see that cancellation is permissible over the whole range of m 3 considered. When the CP 
violating phase a = we see that the leading order term (sin 2^12^23034) vanishes and as a 
result for large values of s§ 4 the cancellation is not possible because the term with coefficient 
A becomes large ((^(lO" 1 )). For non zero values of the CP violating phase a this leading 
order term is non zero and its contribution will be significant. So in this case high values 
of #34 are also allowed because now the leading order and the term with coefficient A will 
be of same magnitude. When we fix s| 4 = 0.06 and a = the region where m 3 is small is 
disallowed (panel (c) blue/dark region) but when a varies within its full range the disallowed 
regions become allowed (panel (c) cyan/light region). When s| 4 approaches its upper limit 
the A term having £ becomes very large and cancellation is not possible even for non zero 
values of a which can be seen from panel (d). However, when a = it/2 very small values 
of s| 4 cannot give cancellation as the leading order term becomes large (panel (e)). s§ 4 has 
to be > 0.01 for the term to vanish which can be seen from panel (f) where we plotted the 
correlation between a and s\ A for \m eT \ = 0. 
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Figure 7: Plots of \m eT \ for inverted hierarchy with lowest mass 777.3. The panel (a) corresponds to 
three generation case. In (b) all the parameters are varied in their full allowed range (3+1). The 
panel (c) , (d) is for specific value of #34 and a with all the other parameters covering their full 
range. The panel (f) shows correlation between a and s| 4 . 



D. The Mass Matrix element m 



/</■< 



This diagonal neutrino mass matrix is given as 
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m m = e 2i ^- s ^c 2 u m 4 s 2 M (29) 
+ e 2 ^ + ^m 3 (c 13 c 2A s 23 - e^-^- s ^s 13 s lA s 24 ) 2 
+ mi{-c 23 c M s 12 + c 12 (-e iS ^c 24 s 13 s 23 - e< 5 ^ 5 ^c 13 s u s M )} 2 
+ e 2ta m 2 {c 12 c 23 c 24 + s 12 (-e iS ^c 24 s 13 s 23 - e^ 5 ^ c X3 s 1A s 2A )} 2 

This expression reduces to its three generation case if the mixing angle 9 24 vanishes. Also 
we can see from the expression that there is no dependence on the mixing angle #34. Using 
the approximation in Eqs. ([2]) this element can be simplified to the form 

\m m \ a I^A^^^v^+e^^^ (30) 

- 2A Cl2 c 23 e i ( 5l3+2a )^ 12S2 3Xi3 

+ X 2 {e 2 ^ +a) ^fCsl 2 s 2 23 x\ 3 + eM^-fe) ^X24 

- 2e 2tQ v / Cci 2 c 23 Si 2 Xi4)X24}}|- 

For a case of vanishing Majorana CP phases and the Dirac phases having the value ir this 
this element vanishes when 

s 23 + C i2 c 23\/C + ci 2 Si 2 sin26» 23V /CAxi3 (31) 

+ ^ (512523^X13 ~ C 23 Sm26 12 ^/CxuX24 + \flxll) = °- 

We know that for the case of 3 generation, the elements in the \x — r block are quite large 
and cannot vanish for normal hierarchy. In panel (a) of Fig. 8 we can see that cannot 
vanish in small ni\ region for s 24 = which is indeed the 3 generation case. This is because 
the magnitude of the first two terms in Eq. f l3T|) is quite large in this case, ~ O (10 _1 ) 
and for cancellation to occur the term with coefficient A 2 has to be of the same order and 
this is not possible when s\ A is small. However when s| 4 varied in its full allowed range the 
contribution of the sterile part is enhanced and this can cancel the active part as can be seen 
from panel (b). Now to understand the dependence of with # 24 , if we increase s 24 from 
its lower bound we can see that the two terms become of the same order. So there will be 
region in the limit of small mi for which this element vanishes (panel (c)). We see in panel 
(d) of Fig. 8 that when # 24 acquires very large values, than the magnitude of the A 2 (v^xD 
term becomes large thus, leading to non cancellation of the terms with the leading order first 
two terms. Hence, the region with very small rrii is not allowed. Using the approximation 
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Figure 8: Plots of vanishing Ulna for normal hierarchy for different values of #24 when all other 
mixing angles are varied in their 3a ranges and CP violating phases are varied in their full range. 



for inverted hierarchy the element m MAt becomes 



\ m nn\ ~ 



^ r77, 23{ C 23( S 12 + C 12 e i 

+ A 2 [sm2^ 12 C2 3 e i(5l4 " 524) (l - 



1 

2' 



-Asin2# 12 sin2#23e i<5l3 (l 



e 2iQ )xuX24 



sl 3 e^(cl 2 + e^sj 2 )xt 3 



e 2ia )Xi3 



(32) 



+ e 



2i(7+5i4-5; 



For a case of vanishing Majorana phases and Dirac phases having value 7r, this element can 
vanish when 



(33) 



In panel (a) of Fig. 9 we plotted |m w | for = to reproduce the case 3 generation whereas 
in panel (b) all the parameters are varied in their allowed range in 3+1 scenario. In both 
the case we can see that cancellation is possible for full range of m 3 . It can be noticed that 
unlike normal hierarchy here cancellation is possible for small values of s\ A it because in this 
case all the terms are of same order and there can always be cancellation. However, if we 
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Figure 9: Plots of vanishing m w for inverted hierarchy with lowest mass rn.3. Panel (a) for 3 
generation case (b) all the parameters are varied in full allowed range (3+1). Panel (c) and (d) 
are for specific values of a and s| 4 are taken with all other parameters covering their full range. 

put a = then the term A (sin 2012S23C23X13) drops out from the equation and the leading 
order term can not be canceled for small values of s\ A . It can be seen from panel (c) that 
for s| 4 = 0.002 and a = the regions where is small is not allowed. As the value of 6*24 
increases there is the possibility of cancellation of terms for all the values of a as can be 
seen from panel (d) where we plot |m MM | with the lowest mass for s\ A = 0.02 when all the 
other mixing angles are varied in 3a range and CP violating phases are varied in full range. 
Now if we keep increasing S24 then A 2 term will become large and the chance of cancellation 
will be less. 
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E. The Mass Matrix element m MT 

The (2,3) element of M v in the flavor basis becomes quite complicated in the presence of 
an extra sterile neutrino. The expression is 

m MT = e^ 2Sl ^ +2 ^c 2 u c 24 m 4 s 24 s u (34) 
+ e 2 ^ + ^m 3 (c 13 c 2 ,s 23 - e l ^-^- s ^s 13 s u s 24 ) 

{— e t<y&1A C 24 Si 3 S\ 4 S 34 + £43(023034 — e %&2i S 23S24S34)} 

+ mi{-C 23 C24Sl2 + Cl2(-e i<5l3 C 2 4S 13 S23 - e 1 ^ 4 ' 5 ^ C 13 S 14 S 2A ) } 
[-Sl 2 (-C 3 4S 2 3 - e i<524 C 23 S 24 S34) 

+ c 12 {-e l5l4 c 13 c 24 s u s 34 - e iSl3 s 13 (c 23 c 34 - e i<524 s 23 s 24 S34)}] 
+ e 2jQ m 2 {c 12C23 c 2 4 + s 12 (-e i5l3 c 24Sl3 s 2 3 - e 1 ^ 5 ^ c 13 s lA s 24 )} 

[Ci 2 (— C34S 2 3 — e i<524 C 2 3S 2 4S34) 
+ •5i 2 { — e 4<5l4 C 13 C 2 4Si4S34 — e^ 13 S13 (c 2 3C34 — e^ 24 S 23 S 24 S 34 )Y\. 

It reduces to the 3 generation case when # 24 = #34 = 0. In the normal hierarchical region 
where m\ can assume very small values and can be neglected, using approximations in Eqs. 
(0 [E2 E} we get 

K T | « | v / A^{c 2 3C34(e 2j(/3+5l3) -e 2 -v / C4)^3 (35) 

- A[c 12 C3 4 e 4{2a+5l3) v^^ 

- ^{e^-^y/lxu - c 12 c 23 e 2iQ ^2x14)334] 

+ A 2 [v^e l(2Q+5l3) (e l5l4 s 12X i4 + 2 Cl2C23 e i<524 x 2 4)si 2S2 3S34Xi3 
+ e^(e^ a - 5 ^c l2 c 34 ^Csi2X2A ~ e^ +5 ^ Xl3 s 3A ) X iAS 2 , 

+ C 2 3C34e 2j ^ 13 ) S 2 2X 2 3V / C^23]}|. 

To see the order of terms we consider the case where Majorana CP phases vanish and Dirac 
phases have the value 7r this element becomes negligible when 

C23C34S 23 (1 - C 2 2 a/C) + -M(ci2C 34 Si2V / CXl3) C0S 2 #23 (36) 
+ X24S 3 4(S 2 3 + c\ 2 C 2 23 \/C,) + 534(^X24 + Ci 2 C 23 Si 2 

VCxu)} 

+ A 2 { s l 2 Xl3 s 23 s 34 12Xl4 + 2c i2 C 2 3X24) + Xl4S 2 3(Ci 2 C 3 4S 2 2 S 23 
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Figure 10: Plots of vanishing m^ T for normal hierarchy (a) for vanishing #34 and #24- In panel (b) 
all parameters are varied in their full allowed range. Panel(c, d) are for specific values of #34 when 
all other mixing angles are varied in their full range. 

Being an element of fir block, m MT shows the same behaviour that of m MAt in normal hierarchy. 
In panel (a) of Fig. 10 we plotted |m Mr | for s| 4 = Sg 4 = to coincide with the 3 generation 
case and we can see that cancellation is not possible in hierarchical region. However, when 
all the parameters are varied in their allowed range in panel (b) it get contribution from the 
sterile part and cancellation is always possible. It can also be seen from panel (c) of Fig. 
10 that s| 4 = is also unable to give cancellation in the region when mi is small and the 
figure is quite similar to that of 3 generation case. However, as this active sterile mixing 
angle becomes larger there is always a possibility of allowed region towards the lower values 
of mi as is evident from panel (d). This is because for the vanishing value of 6*34 the terms 
with A and A 2 become very small and cannot cancel the leading term O (10 _1 ). It can also 
seen that in this case (i.e s| 4 = 0), there is no X24 term in Eq. (3.31) and this is why the 
figure is somewhat similar to the 3 generation case. However, when # 34 increases these two 
contributions become large and thus cancellation becomes possible. For the case of inverted 
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hierarchy where can have very small values, m MT becomes 

\m^ T \ » |y / Am^{-C23C 3 4S23(ci 2 e 2jQ + s? 2 ) (37) 
+ A[c 12 s 12 (l - e 2ia )(c34Cos2^ 3 e i5l3 xi3 + c 23 s u e iSl * X u) 
+ s, A {e^^^l- cl^{s\ 2 + c 2 2 e 2 -)} X24 ] 
+ A 2 [ C2 3C34S 23 e Ml3 (c 2 2 + e 2ja S 2 2 )x 2 3 

+ c 12 s 12 s 23 (e 2ia - l)(c 3 4e l(,5l4 - 524) xi4 + 2 S3 4C 23 e 4{5l3+524) Xi3)X24]}|. 

To get an idea about the magnitude of the terms we take the vanishing Majorana phases 
and Dirac CP phases to be of the order tt. The expression in this case for vanishing m^ T 
becomes 

-C23C34S23 + A(s 3 4X24(c 23 - v 7 ^)) - A 2 (-s 34 Xi3Xi4 - C23C34X13) = (38) 

In panel (a) of Fig. 11 where Im^l for 3 generation is plotted we can see that unlike 
m w there is no cancellation in small m 3 region but when plotted for the full range it gets 
contribution from the sterile part and there is cancellation for the full range of (panel 
(b)). Clearly the cancellation of the terms do not become possible for small values of #34 in 
strict hierarchical region. This case is similar to the three generation case in IH (cyan/light 
region, panel (c)). This is because for s| 4 = the contribution of s 24 comes from the A 2 
term. If we put the CP violating phase a as zero then cancellation is not possible for whole 
range of m 3 (blue/dark region panel (c)). However, as the value of s 2 4 increases all the terms 
in the above equation becomes of the same order and cancellation for very small values of 
m 3 is possible (panel (d)). 

F. The Mass Matrix element m TT 

This element is related to m w by the fi — t symmetry. As discussed earlier, in the limit 
when # 2 4 and #34 are not very large, the two mixing angles #34 and #24 will behave same in 
the textures related by fi — r symmetry. The (3,3) element of the neutrino mass matrix in 
the presence of one sterile neutrino is given as 
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Figure 11: Plots of vanishing m^ T for inverted hierarchy (a) for vanishing #34 and #24 (3 
generation). In panel (b) all parameters are varied in their full allowed range (3+1). Panel (c) and 
(d) are for specific values of #34 and a when all other mixing angles are varied in their full range. 



Tfl TT — 6 ' C^C 2 ^Tn^S 

D 2i(<5 13 +m r J(S 14 -S 13 ) 



(39) 



+ e^ +IJ} m 3 {e^~ d ^c 2i s 13 suS34 + c 13 (c 23 c u - e*« s 23 s u s u )Y 

+ "ii[-si 2 (-c 3 4S 2 3 - e l524 c 23 s M s 34 ) 

+ Ci 2 {— e lSl4 Ci 3 c 2 4Si4S 3 4 — e iSvi Si 3 (c 23 c 3 4 — e l<524 S23S24S34) }] 2 

+ e 2ja m 2 [ci2(-C34S23 - e l&2i c 23 s 2i s 3i ) 

+ s 12{ — e t(Sl4 Ci3C24Si4S34 — e** 13 S\ 3 (C23C34 — e^ 24 ^^^)}] 2 - 

It reduces to the 3 generation case for # 34 = 0. Using the approximation for normal hierarchy 
in Eqs. P US E]) this becomes 
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\m TT \ « | y / Am^{c 23 C34S23(e 2i/3+<513 - c\ 2 ^e 2ia ) + A{-e*( 2a +^ (40) 

Sl 2 C 12 C 3 4 COS 2^3X13 - ^^S^ {s l2X l^ + C 12 C 2zX 2^ ) + S 3 4X24(-^ 3 e 2 ^ /3+<5l3 ) +i524 

- + A 2 {v^s 2 2S23 Xi3e l(2a+5l3) (c 2 3C34Xi3e^ 3 + s 3 4Xi4e i514 ) 
+ v / Cci 2 si 2 S23X24e 2ia (2c 23 S34Xi3e i{5l3+524) + c 34 Xi4e i(5l4 ^ 24) ) 

- ^23 S 34Xl3Xl4e i(2/3+5l3+5l4) }}|. 

To get an idea of the order of the terms we consider the vanishing Majorana phases and the 
Dirac phases having the value equal to ir. This element vanishes when 

C 23 C 3 4S23(1 - C? 2 a/C) + A{a/CS12Ci 2 C 34 COS 20 23 Xl 3 + a/Cc 12 C 23 S3 4 (Si 2 Xi4 (41) 
+ Ci 2 C 23 X 24 ) - S 34 X24(S 23 - V 7 !)} + A 2 {v / Cs 2 2 S 23 X13(c 2 3C 3 4X13 + S34X14) 
+ \/C C i2 s 12S 2 3X24(2c 2 3S34Xl3 + C 3 4Xu) + S 23 S 3 4Xl3Xl4} = 0. 

It is found that for vanishing 6*34 which is the case for 3 generation is disallowed for small 
mi as can be seen from panel (a) of Fig. 12. This is the generic behaviour of a element 
belonging to the \x — r block in normal hierarchy which we mentioned previously. This is 
because for 634 equal to zero the leading order term is large (O (10 -1 )). Here the term with 
A 2 is quite small (10~ 3 -10 -4 ) and hence will not have very significant role to play. Thus, only 
term with coefficient A can cancel the leading order term. However, for vanishing 634 this 
term is small O (10~ 3 ), thus not leading to any cancellation of the leading order term. In 
panel (b) when all the parameters are varied in their 3a range we can see that cancellation 
is possible over the whole range of mi (3+1 case). Now, when #34 starts increasing from 
its lowest value there exist a region for intermediate values where both the terms become 
approximately of the same order and hence there can be cancellations (panel (c)). Towards 
very large values of 634 the term with coefficient A becomes larger than the leading order 
term due to which this element cannot vanish. For the cancellation very large values of mi 
is required as can be seen from panel (d) of Fig 12. For the case of inverted hierarchy where 
m 3 approaches small values we get the expression 

m rr « c 2 4S 2 3 ( C 2 2 e 2 - + s\ 2 ) + e 2 ^) ^34 (42) 
+ 2A[(e 2ia - l)ci 2 c 34 Si 2 s 2 3(c 2 3C34Si 2 s 2 3(c 2 3C 34 e 2i5l3 Xi3 + e 2 '' 5l4 s 34 Xi4) 
+ 2c 23 C34e 2j524 S2 3S34(c 2 2 e 24a + S 2 2 )X24] 

+ A 2 [(c 2 2 + e 2M S 2 2 ){ C2 3C3 4 Xi3^ 13 (c 2 3C34Xi3e 4 ' 513 + 2xi4S34e* 514 ) + e 2i<5l4 X 2 4 4) 
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Figure 12: Plots of vanishing m TT for normal hierarchy with lowest mass mi. Panel (a) 
correspond to three generation case. In panel (b) all the mixing angles are varied in their full 
allowed range of parameters (3+1). Panel (c) and (d) are for some specific values of #34. 



+ (c 12 e + s 12 )c 23 e 24 X24 S 34 

+ 2s 12 {e 2ta - l)e i<524 (c34Xl3COs2^3e^ + C^Xl^SsiX^]- 



For vanishing Major ana CP phases and Dirac phases having the value equal to it this ex- 
pression becomes 



m TT « -C23C34S23 + AS 34 X24(C23 - V^0 + A 2 S23Xl3(c 2 3C 3 4Xl3 + S 3i Xu) (43) 



In panel (a) of Fig. 13 we reproduced the 3 generation behaviour by plotting \m TT \ for 
s| 4 = and in panel (b) all the parameters are varied randomly (3+1). In both the cases 
we can see that cancellations are possible for the whole range of m 3 . For s| 4 = all the 
terms are of same order and cancellations are always possible. But if we put a = then one 
term with coefficient A and another term with coefficient A 2 drops out from the equation 
and then small values of s| 4 can not cancel the leading order term any more. This can be 
seen from panel (c) where cancellation is not possible for lower 777,3 region. However when 
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Figure 13: Plots of vanishing m TT for inverted hierarchy with lowest mass m^. Panel (a) 
correspond to three generation case. In panel (b) all the mixing angles are varied in their full 
allowed range of parameters (3+1). Panel (c) and (d) are for some specific values of #34 and a. 

s| 4 increase to a value of about 0.02 this element can vanish (panel (d) the cyan region). We 
see that when #34 increase towards its upper bound the A term becomes large 0(1). Hence, 
the other terms are not able to cancel this term so we do not get small m 3 region allowed 
(Panel (d), blue region). 



G. The Mass Matrix elements m es , m„ s , m TS and m ss 

The elements m es , m^, m TS and m ss are present in the fourth row and fourth column in 
the neutrino mass matrix. They are the new elements that arises in 3+1 scenario due to the 
addition of one light sterile neutrino. The expressions for m es and m MS are given by 
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m es = e i{2l+5li] ^014024034^4514 (44) 

+ e i(2/3+5l3) c 1 4m 3 s 1 3{-e 4(5l4 - 5l3) c 2 4C34Si3S 1 4 + c 13 (-e t5l *c u s 23 s 2 ± - c 23 s u )} 

+ Ci2Cl3Ci^mi[—Si2( — e lS24 C2 3 C 3 4S24 + S23 s 34) 

+ Ci2{ — e 4<5l4 C 13 C24C34Si4 — e^ 13 S13 (— e^ 24 C34S23S24 — C23S34)}] 

+ e 2ia Ci 3 Ci 4 m2Si2[Ci2(-e i<524 C23C34S24 + S23S34) 

+ Si 2 {-e J<5l4 Ci3C24C 3 4Sl4 - e i5r3 S 13 {-e l&2i C 3i S 23 S 2i ~ C23S34)}]. 



m„ s = e^+^c^cuCum^SM (45) 
+ e^ +5 ^m 3 (c 13 c 2i S2 3 - e^- s ^) Sl3 s lA S24) 

{-e i(5l4 - 5l3) C24C 3 4Si 3 si4 + c 13 (-e lS24 c 34 s 23 s 2 4 - C23S34)} 

+ m 1 {-C23C24S 12 + C 12 (-e i5l3 C24Sl3S 2 3 - e^ 14 "* 24 Wl4S 24 ) } 

[-Si 2 (-e^ 24 C23C34S 2 4 + S23S34) 
+ Cn{ — e l<5l4 Ci3C24C34Si4 — e i5l3 S 13 ( — e l<524 C34S23S24 — C23S34)}] 

+ e 2ia m2{c 12 C23C24 + s 12 (-e i5l3 C24S 13 S23 - e j(5l4 - 524 c 13 s 1 4S24)} 

[Cl2(-e" S24 C23C 34 S24 + S 23S 34 ) 
+ — e i<5l4 Ci3C24C 3 4Si4 — € %&13 S\ 3 (— e^ 24 C34 S 23 s 24 _ c 23 s 34)}]- 

Though the equations seems very complex, one can easily understand the properties of these 
elements by just looking at the m 4 terms. The m 4 term in m es is proportional to si 4 . So in 
general it is quite large 0(1). For this element to become negligible very small values of s 2 4 
is required. But as this angle is bounded by the SBL experiments, complete cancellations 
never occurs for both normal and inverted hierarchy (Panel (a), (b) of Fig. 14). Similar 
predictions are obtained for where element cannot vanish as there s| 4 has to be negligible 
which is not allowed by the data as can be seen from Panel (c), (d) of Fig. 14. 
For the element m TS the scenario is quite different. 
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Figure 14: Correlation plots for vanishing \rrik s \ for both normal and inverted hierarchy. In these 

plots all the mixing angles are varied in their 3 a allowed range and the CP violating phases are 

varied in their full allowed range. 
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m TS = c? 4 c4c34e 2i(5l4+7 W 3 4 (46) 
+ e 2 ^ +s ^m 3 {-c 24 c 34 e^- s ^s 13 s u + c 13 (-c 23 s u - c u e l5 * 4 s 23 s 24 )} 

{ — C^e 1 ^ 14 Sl3 ' S13S14S34 + C\ 3 (c 23 C 3 4 — e t52i s^s^s^)} 
+ mi[-Sl2(-C23C 3 4e l<524 S 2 4 + S23S34) 

+ c 12 {-c 13 c 24 c u e iSl4 s 14 - e i5is s 13 (-c 23 s 34 - c 34 e i<524 S23S34)}] 

[-Si 2 (-C 34 S23 - C23e l<524 S 34 S24) 

+ c 12 {-c 13 c 24 e i5l4 s u su - e i5ls s 13 (c 23 c u - e^ 24 ^^^)}] 
+ e 2iQ m 2 [ci 2 (-C23C34e^ 24 S34 + S23S34) 

+ S\ 2 {— Ci 3 C 2 4C 3 4e lSl4 Si4 — e i<5l3 Si3(— C23S34 ~ c 34 e * <524 '523 s 24)}] 

[C1 2 (-C 3 4S23 - C23e J<524 S 34 S24) 
+ s 12{~ Cl3 c 24e l<5l4 Si4S34 — e^ 13 S 13 (C23C34 — e^ 24 S23S34S14) }] . 

In this case the m 4 term is proportional to #34 and there is no lower bound of it from the 
SBL experiments i.e. it can approach smaller values. As a result the term with m 4 can be 
very small. Thus this matrix element can possibly vanish in both hierarchies for whole range 
of the lowest mass (Panel (e), (f)). 

The (4,4) element of the neutrino mass matrix is given as 

m ss = cf 4 44e 2! h +fe 'm 4 (47) 

+ e 2i(/3+<5l3) m 3 {-C24C34e i(5l4 - <5l3) s 13S i4 + c 13 (-C2 3 S34 - c 34 e J<524 S 23S24)} 2 

+ mi[-si2(-C23C 3 4e J<524 S24 + S23S34) 

+ Ci2"{— c i3 c 24C34e J<5l4 Si4 — e l8u S13 (— C23S34 — c 3 4e tS24 S23S24) }] 2 

+ e 2ia m 2 [c 12 (-c 23 c 34 e i524 s 24 + s 23 s 34 ) 

+ S\ 2 { — Ci3C24C34e* <5l4 S i 4 — e lSl3 Si 3 ( — C23S34 — C34e l<524 S23S34)}] 2 

The WI4 term for m ss is proportional to c^c^c 2 ^. One can see that this term is of order 
one as a result this element can never vanish as is evident from panel (g, h). 



IV. CONCLUSIONS 



In this paper we analyze systematically the one-zero textures of the 4x4 mass matrix 
in presence of a sterile neutrino. Assuming neutrinos to be Majorana particles, this is a 
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symmetric matrix with 10 independent entries. We analyze if the current constraints on 
oscillation parameters allow each of these to assume a vanishing value. We also study 
the implications and correlations among the parameters when each matrix element is zero. 
We expand the mass matrix element in terms of a parameter A with suitable coefficients 
Xi4;X24 and X34 corresponding to the mixing angles 8 13 , 8 M and # 2 4- This is motivated 
by the observation that these angles are of order ~ 0.02. These expressions facilitate the 
analytic understanding of the numerical results presented in the different plots. We study 
the vanishing condition as a function of the lowest mass m\ (NH) or m 3 (IH) by varying the 
lightest mass in the range 0.0001 - 0.5 eV. 

We find that \m ee \ = is possible for NH only for higher values of the smallest mass 
mi while for IH it is possible even for lower values. This is in sharp contrast with the 3 
generation case where complete cancellation can never take place for IH. 

|m e/ J can vanish over the whole range of the smallest mass for both 3 and 4 neutrino 
scenarios. However for larger values of the mixing angle s\ A cancellation is not achieved for 
smaller m\ for NH. For IH the the cancellation condition depend on the Majorana phase a 
and the mixing angle #24- We obtain the correlations between these two parameters required 
for making this element vanishingly small. 

Cancellation is achieved for the element m er for the full range of the lowest mass in the 
3+1 scenario. The element m eil is related to the element m eT by fi — r symmetry. However 
unlike three generation case 823 in these textures are not related simply by 823 = (tt/2 — 823). 
The mixing angles 824 and 834 are also different in these two textures in general. However 
for small values of 6*24 we get #24 = 834 in these textures. Consequently the role played by 
#24 for m eAt is played by 834 in this limit. Thus in this case cancellation is not achieved for 
larger values of s| 4 in the hierarchical regime for NH. For IH we obtain correlations between 
a and sin 2 834 for fulfilling the condition for cancellations. 

The elements and related by \x — t symmetry. For these cases, cancellation 

is not possible in the hierarchical zone for IH in the 3 generation case. However the extra 
contribution coming from the sterile part helps in achieving cancellation in this region. For 
IH one can obtain correlations between the Majorana phase a and the mixing angle 824(834) 
for |m MM | = 0(|m rr | = 0). 

For m Mr element cancellation was possible for three generation case only for higher values 
of the lightest mass. However if one includes the sterile neutrino then this element can 



31 



vanish over the whole range of the lightest neutrino mass considered. 

With the current constraints on sterile parameters it is not possible to obtain m ss = 
while m es and can only vanish in the QD regime of the active neutrinos. On the other 
hand m TS can assume vanishingly small values for the lowest mass in the range 0.0001- 0.5 
eV, since the angle #34 has only an upper bound from current data and hence can be very 
small. 

The above results can be useful for building models for light sterile neutrinos and shed 
light on the underlying new physics if future experiments and analyses reconfirm the expla- 
nation of the present anomalies in terms of sterile neutrinos. 
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